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Abstract 

A theory is offered for the balhstic 4-terminal Josephson junction. 
The studied system consists of a mesoscopic two-dimensional normal 
rectangular layer which is attached in each side to the bulk supercon- 
ducting banks (terminals). The relation between the currents through 
the different terminals, which is valid for arbitrary temperatures and 
junction sizes, is obtained. The nonlocal coupling of the supercurrents 
leads to a new effect, specific for the mesoscopic weak link between two 
superconducting rings; an applied magnetic flux through one of the 
rings produces a magnetic flux in the other ring even in the absence 
of an external flux through the other one. The phase dependent dis- 
tributions of the local density of Andreev states, of the supercurrents 
and of the induced order parameter are obtained. The "interference 
pattern" for the anomalous average inside the two dimensional region 
can be regulated by the applied magnetic fluxes or the transport cur- 
rents. For some values of the phase differences between the terminals, 
the current vortex state and the two dimensional phase slip center are 
appeared. 



1 



1 Introduction 

The Josephson multiterminal junction presents the microstructure in which 
the weak couphng takes place between several massive superconducting banks 
(terminals) [1-3] . Compared with the conventional (2-terminal) Josephson 
junctions such systems have additional degrees of freedom and correspond- 
ing set of the control parameters. With the result, for example, the current- 
or voltage- biased and the magnetic flux-driven regimes can be combined 
in one multiterminal junction. The specific multichannel interference effects 
were studied theoretically and experimentally in the novel superconducting 
device, the 4-terminal SQUID controlled by the transport current [5-7]. Re- 
cently another system based on a Josephson 4-terminal junction was studied 
[|]. It consists of two superconducting rings, each interrupted by a Josephson 
junction, which are at the same time weakly coupled with each other. The 
macroscopic quantum states of such composite system can be regulated by 
the difference of the magnetic fluxes applied through the rings, in analogy 
with the phase difference between two weakly coupled bulk superconductors. 
The nonlinear coupling via the Josephson 4-terminal leads to the cooperative 
behaviour of the rings in some region of the applied magnetic fluxes, which 
was called |^ magnetic flux locking. 

The 4-terminal junction, which was studied in Refs. [5-8], is a system 
of short microbridges going from a weak point to massive superconducting 
banks (Fig. la). The order parameter (both its amplitude and phase) in the 
common centre is a function of the currents through all the microbridges. 
The supercurrent flowing into the ith bank is determined by the phases of 
the order parameter yjj, (i = 1, ..4) in all the banks f^: 



The relation (1) was obtained in the frame of the Ginzburg-Landau approach, 
which is valid for temperatures T close to the critical temperature T^. As 
it was pointed in the Ref. [^j, the macroscopic interference effects due to 
coupling of supercurrents in different terminals are not restricted by the 
special kind of the 4-terminal junction (Fig. la). In fact, any mesoscopic 4- 
terminal weak link will produce a coupling similar to the relation (1). In the 
present paper, the microscopic theory of the mesoscopic ballistic 4-terminal 




(1) 
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junction is developed. We consider a Josephson weak coupling through the 
two-dimensional normal layer which is connected with four bulk supercon- 
ducting terminals as it is shown in Fig. lb. Such a S-2DEG-S structure was 
experimentally realized in Ref. for the case of two terminals. It was 
shown in that this new class of fully phase coherent Josephson junctions 
demonstrate the nonlocal phase dependence of mesoscopic super cur rents. We 
study the coherent current states in such a 4-terminal structure within the 
quasiclassical equations for transport-like Green's functions. The relation 
between the currents in the different terminals, which is valid for arbitrary 
temperatures and junction sizes, is obtained. The structure of current car- 
rying states inside the mesoscopic 4-terminal junction presents the interest 
itself. As it is well known (see e.g. |Tn|), in ballistic Josephson junction with 



direct conductivity the supercurrent flows through the local Andreev lev- 
els. In multiterminal case considered here, the spatial distribution of current 
density and of the order parameter, and hence the phase-dependent Andreev 
levels, are determined by the phase differences between all terminals. Thus, 
they can be regulated by the external control parameters, i.e. the transport 
currents and (or) the applied magnetic fluxes. In Section 2, we present the 
description of the system and formulate basic equations and boundary con- 
ditions. In Section 3 the current-phase relations analogous to (1) are derived 
for the cases of small (as compared to the coherence length) and also arbi- 
trary junction sizes. The spatial distributions of the supercurrent density 
and of the induced order parameter are studied in Section 4. 



2 Model and basic equations. 

The studied system consists of 4 bulk superconducting banks which are con- 
tacted with 4 sides of rectangular two-dimensional (2D) normal layer having 
the length L and width W (see Figs, lb and 2). The sizes L and W are 
supposed to be much larger than Fermi wave-length \p = h/pp ■ To study 
the stationary coherent current states in 4-terminal ballistic junction we use 



the Eilenberger equations |11| for .^-integrated Green's functions: 



vp|-G+[cuf3 + A,G] = 0, (2) 
or 



where 
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is the matrix Green's function, which depends on the Matsubara frequency 
u, the electron velocity on the Fermi surface and the coordinate r; 




is the superconducting pair potential. For the self-consistent off-diagonal 
potential A(r) and current density j(r) we have the expressions 

A(r) = A27rT ^ (^) , (3) 

a;>0 

j(r) = -47rzeiV(0)T^(v^(7..). (4) 

a;>0 

They determine the induced order parameter \E'=A/A and the 2D current 
density in the normal layer; A^(0) = (...) is the averaging over directions 
of 2D vector vp, A is the constant of electron-phonon coupling. 

The equations (2) are supplemented by the values of A and Green's func- 
tions in bulk banks far from the S-N interfaces: 



Aj = Ao('ricosv3i - f2sin(^i), = -^^^^=i^, i = l,.A. (5) 

+ A§ 



We solve Eqs. (2) by integrating over the "transit" trajectories of the ballistic 



flight of electrons from one bank to another |jT2|. These trajectories (charac- 
teristics of the differential Eqs. (2) ) are straight lines along the direction of 
electron velocity (see Fig. 2). In the bulk superconducting banks the order 
parameter can be taken as the constant value (5) up to the S-N interface. 
In contrast to the case of 2D banks these "rigid" conditions for A |Q |12| 
are valid for arbitrary sizes L and W compared with the coherence length 
^0 ~ vp/Aq, and not only for L,W C,^. At the same time, the Green's 
function along the given transit trajectory varies in a distance of about 
where approaching the S-N interface. 
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Let us introduce the time of flight along the trajectory, vp-^^ = ^, 
ti < t < oo where t = ti corresponds to the point on ith S-N boundary 
and t = oo to the point inside ith bank far from the S-N boundary. Then 
the general solution of Eqs.(2) inside the ith bank satisfying the boundary 
conditions (5) will be 



+ Ci[Ao'r3 — {ujcosifi + zsz(7n(vi7nj)fi sin (y^j)-?-! (6) 
+ sin(y9j — isigniy F'!a.i)VL cosv9j)'r2]e^^^^*~*'''. 
Here rij is the outer normal to the ith side of the rectangular boundary and 



VL = ^u'^ + Aq. The arbitrary constants Ci must be found by matching 
of Green's functions at in-coming and out-going points at S-N boundaries 
with the solution inside the normal layer along the trajectory which connects 
these points (see Fig. 2). We consider here the simplest case when only 



Andreev reflection 13 occurs at S-N interface. In more realistic case, when 



usual reflection (e.g. due to the potential barrier) or interface roughness are 



present, more general matching conditions must be used (see W4 



3 Current-Phase relations 

Inside the normal layer (A = 0), the Eilenberger equations can be solved 
analytically. If we classify the electronic trajectories inside the normal layer 
according to the sides at which they come in and go out, then the solution 
of eq. (2) can be written as 

Gi^jit) = ^fg + ^{cosh [2iu{t - ti) - iipi]ri - i sinh [2u{t - U) - i(pi]f2} 

-|-y4j_^j{AoT3 — [uj cosh {2uj{t — ti) — iifi) + Q sinh {2uj{t — ti) — iLpi)]fi 

+ i[Q cosh {2uj(t — ti) — iifi) + uj sinh (2ijj{t — ti) — iipi)]f2}, (7) 

where Gi^j(t) is the matrix Green's function along the trajectory originat- 
ing in the ith side and extending to the jth side (see Fig. 2). We denote 
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this trajectory by i j- Matching (7) with solution in the banks (6), the 
corresponding Ai^j is obtained 

uj sinh {ujtji + i^) + O cosh {ujtji + i^) ' 

where tji — tj — ti and ipji = ipj — ipi. Prom (7) and (8) we have the expression 
for matrix Green's function Goj{p, Vp) as a function of the coordinate p 
(E is the region of 2D rectangular weak hnk) and the direction of Vp- In fact 
we can write 

G^{p,^f) = Gi^j, for vp e 9ij{p}, (9) 

where we have introduced as the angle in which all i ^ j trajectories, 

passing through the point p, are confined (see Fig. 2, eq. (Al) in Appendix). 
The diagonal and off-diagonal terms of Gi-^j have the forms 



LO cosh i^ij, + /^) + 9. sink {^ij, + /^) 



(10) 



n cosh {ujtji + i^) + a; sinh {totji + i^) 

fi^j " + ~ A^i] [-2a; (^ - k) + m] 

^ Aoexp(a;t,-, + i^) g-2a,(t-ti) /^^x 

Q cosh (a;iji + i^) + a; sinh (a;iji + i^) 

In the hmit L, 1^ << the expressions (10) and (11) for Green's functions 
are simplified and we have 

_ a;Q + i|Ao^sin((^j^) 
a;2 + Ao'cos2(^) ■ 

f _ Ao 

Qcos(^)+^cc.sin(^)^''P^ 2 ^" ^ ""^ 

We can obtain the retarded and advanced Green's functions, G'^'^(e), 
by analytical continuation of Matsubara Green's function G{uj) (eqs. (9)- 
(13)). The poles of diagonal component of the retarded Green's function. 
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g^{e, PjVf), determine the energies of local Andreev states in the system. 
The local density of states in the normal layer is given by the furmula 

Af{e, p) = N{0) < Re g{uj = -ie, p, v^) > . (14) 

Using the expressions (9), (12) and the fact that, Oij{p) — Oji{p), in the case 
of small junction, we obtain: 

^(e, P, Wi}) = N{0) Y,<Reg{'^ = -^e, P, vf) >0ij 

= iV(0) J2 %(p)^e g.^jiou = -ie) 

= N{^)^^ij{p)R49i-^ji^ = -^e) + 9j-^i{^ = -i^)] 

i<j 

= 7rAo7V(0)$:%(p) I sin ^ I S{\ e \ -AqCOS^). (15) 

We can also use eqs. (9) and (12) to obtain < \pg > at a point of the 
ith side pi. Then, the resulting expression can be replaced in eq. (4) to 
find the current density j(pi). The calculation of the current density at the 
arbitrary point of the normal rectangular will come in the next section and 
the appendix. Here we calculate the total current /j flowing into the ith 
bank. 

Let us start with the case of a small junction {L,W « ^q). In order to find 
li, we have to calculate the integral Jj = /(^.-j j(pi).dsi, where the integral is 
taken over the ith side of rectangular. 

After calculation of j(pi) from (A3) and (A5) and taking the integral over 
dSj, we obtain for the current /j 

eppAodJ^ . .(Pi-(pj.. .Aocos(^^^) 
li = ^ ^'"^ ^ — 2 ^ ^ 2r ^ ^ ^ 

where d = \/LP~+~W^ and = 7ji, 
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1.1 + k , , , 

712 = 714 = 723 = 734 = 2^ ^i + k^ ~ ^-^^^ 

are geometrical form factors which depend on the width to length ratio k — 
W/L. The positive sign of Jj corresponds to the direction of the current from 
the normal layer to the ith bank. Note that Y^j=i li = 0. 
The formula (16) for current-phase relations generalizes the expression (1) 
to the case of a small mesoscopic 4-terminal junction. It follows from (17) 
that the form factor ^y^j can not be factorized, i.e., presented in the form 
lij = 7j7j) in contrast to the case of relation (1) where jij = This 
essential feature of the current-phase relations reflects the nonlocal nature of 
the supercurrents in the mesoscopic multi-terminal Josephson junction. 

The current-phase relations (16) are valid for arbitrary temperature T. 
In the hmiting cases of T = and temperature close to Tc the expression 
(16) takes the forms 

5 Z^7iiSin( —^) forT^O, (18) 

^ i=i ^ 

/. = ^^^E7.,sin((^,-(^,) fovT^T,. (19) 

In the case of arbitrary lengths L, W we restrict the consideration for the 
temperature close to T^. In this case the current-phase relations similar to 
the expression (19) can be obtained. The difference is in geometrical form 
factors. In fact we have the result 

h = E UK L, W) sin (v., - VP,) (20) 

where the generalized form factors are given by 



741 = 721 = 723 = 743 



, : dv dd cos 

TT^yiTFi-l iarctan(|+j/) (2n + 1)2 
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~ ~ 4 n , .arctand-,) ^ ^ exp [-^(2n + l)] 

742 = 724 = , / ay / dUco&ti > r- , 

(21) 

7l2 = 7l4 = 732 = 734 = 74l(-^ ^ H^, ^ iv, ^ 

713 = 731 = %2{L ^ W,W ^ L,k ^ 
Here ^jv = In the limit L,W « ^q, 'jij reduce to jij. 



4 Spatial distribution of supercurrents and 
induced order parameter 

In this section we will obtain the supercurrent density and the induced order 
parameter at an arbitrary point of the normal layer in the case of small 
junction. At the given point of the normal layer p — xi + 



< vf^ >= ^(< vf^ >0,,. + < vpg >eji) 



i>j 

2 , 



o;^ + Ao COS"' ^ 



Ao sini^j 

i>j - +Ao^COS' 2 

where we have used < VPS' >eij^< vp >eij Qi-^j, < vp >e,i= - < vp >eij 
and fifj^i = g*^j. 

Replacing (22) in the eq. (4) the current density is obtained 

j(p) = 27reiV(0) Ao E < vf >e., sin ^ tanh ( °^^'~ ) (23) 

i>j 

The expression (23) describes the spatial distribution of the current density 
inside the normal layer. In order to find the explicit expression for the coef- 
ficients < Vp >eij in eq. (23), we have to consider 4 different regions in the 
normal rectangular and obtain j(p) in each region separately (see Appendix). 
This calculation has been done in the appendix and the result for j(p) is given 
by (A3) and (A5). Here we write eq. (23) in the more transparent form. Let 
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us introduce 0ij{p) as the unit vector in the direction of the i ^ j trajectory 
passing through the bisector of Oij{p); then, < vj? >0^. can be wrirren as 

Combining eqs. (23) and (24), we obtain 



eppAo^ . Oij^ . (fj-(fii Aocos(^^i2^) 
J(P) = Z^sm {^)9ij{p) sm (^^^ ) tanh [ ]. (25) 



i<j 



In the way similar to what we have done for j(p), the distribution of the 

induced order parameter can be obtained. In this cas we need to calculate 
the average of off-diagonal element of the matrix Green's function, fuii^F, p), 
on the direction of vp: 

< / >= E < / >e..= E (26) 

Replacing (26) in (3) and after the calculation, we obtain for ^(x, y) = 
A(a:,|/)/A: 

^{x,y)^^Y.Q,{x,y)e^^K (27) 

Here ^j(x, y) is the angle by which ith side is seen from the point p = {x, y). 
The angles 9i{x,yys are given by the relations 

0^=Ti- a{x, y) - tt(-x, y), 6*2 = a{-x, -y) + a{-x, y), 
63 = 11 - a{x, -y) - a{-x, -y), 64 = a{x, y) + a{x, -y). (28) 
where the angle 

a{x, y) = arctan (jy^^)' (29) 

is a function of the coordinate (normalized by L) and is shown in Fig. 6. 
Eq.(27) expresses the fact that, inside the ballistic normal layer region the 
linear superposition of four macroscopic wave functions (pair potentials) of 
the banks occurs, where the weight of wave function of the ith bank is de- 
terminened by the geometrical factor 6i{x,y). 
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5 Conclusions 



The present study considers a 4-terminal microstructure based on a new 
class of mesoscopic Josephson junctions [9] which arc fully phase coherent and 
have comparable width and length. The microscopic theory of the stationary 
coherent current states in ballistic multiterminals is developed. 

We have calculated the current- phase relations (CPR), i.e. the total cur- 
rents in each terminal as functions of the phases of the superconducting order 
parameter in all the banks. These relations describe the behaviour of the sys- 
tem influenced by the external transport currents or the applied magnetic 
fluxes. The essential difference between the CPR for mesoscopic (expression 
(19)) and conventional (relation (1)) 4-terminals consists in the structure of 
the coefficients of coupling 7^-. In the mesoscopic case considered here these 
coefficients can not be factorized (presented in the form jij = •ji'-yj) for all in- 
dexes i,j and arbitrary value of the width to length ratio k = W/L. Here we 
only outline the new effect, specific for the mesoscopic 4-terminal junction, 
which follows from such nonlocal coupling of the currents. Let us consider 
the configuration shown in Fig. 3. By using the CPR (19) with given by 
(17), it can be shown that an applied magnetic flux through one of the rings 
produces magnetic flux in the other ring even in the absence of an external 
flux through the other one. The detailed theory of this effect will be reported 
in a separate publication. 

The physical properties of the interior of the mesoscopic 4-terminal junc- 
tion present the interest themselves. The above calculated local density of 
Andreev states, the current density and the order parameter distributions 
depend on the phase differences between the four terminals and can be reg- 
ulated by the applied magnetic fluxes. In particular, for some values of the 
phases </?, 9 and x (see Fig. 3) the "vortex state" inside the mesoscopic 2D 
weak link exists. Figures 4 and 5 present the plots for distributions of the 
absolute value of the induced order parameter and the supercurrent density 
in the case 9 = 7t/2, ip = 37t/2, x = 0. The studying of the structure of 
induced order parameter and local density of states, as well as the dynamical 
behaviour of the system will be the object of further investigation. 
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Appendix 



In this appendix we present expressions for the angles 6ij and the vectors 
< >e, ■ Using the expressions given here, one can calculate the density 
of states A/" and the current density j [see eqs. (15) and (23)]. 

According to the classification of the trajectories in term of origin and 
destination sides, there are 12 different types of trajectories which are 1 — >■ 
2, 1 — > 3, 1 — > 4, 2 — > 3, 2 — > 4, 3 — > 4 and the corresponding reverse of these 
trajectories. For a given point, depending on the position, some of these 
trajectories do not take place. In this respect we can consider four different 
regions in the normal rectangular: / where y < 0,\ y \> k \ x \,{2 ^ 3,3 ^ 4 
and their reversed are absent) 

// where x>0,|y|<A;x,(l— >4,3— >4 and their reversed are absent), 

III where y > 0,y > k \ x |,( 1— >2,1— >4 and their reversed are absent) 
and 

IV where X < 0,\ y \< kx,{ 1 ^ 2, 2 ^ 3 and their reversed arc absent). 

At the given point p, for the absent trajectories we have 9ij = 0, and conse- 
quently the corresponding term in the expressions of J\f and j [eqs. (15) and 
(23)] will vanish. We will calculate j in the given point of the region // and 
then introduce the exchange rules of arguments to obtain it in other regions. 
Consider a point in region //; the possible (non-vanishing) 6ij arc drawn in 
Fig. 6 and, can be expressed in terms 6'i's (given by (28) and(29)) as 

012 — -{Oi + 02 — O3 — ^4) 

^13 — -^{01 — 02 + 03 + 04) 

023 — 2^~^^ + 02 + 03 — 04:) 

024 = 04 (Al) 

Also we can use the relation, < vp >eij= I(g^^) |f^^F(cos^i-|-sin0j), to obtain 

< Vf >ei2 {x,y) = ^{[sm[a{-x,-y)] - sm[a{x , y)]]i 
+[cos[a{x, y)] - cos[a{-x, -y)]]i} 
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< Vf >ei3 {x,y) = —{[sm[a{~x,y)] - sm[a{-x, -y)]]i 
Zn 

+[cos[a{-x,y)] + cos[Q;(-a;, -y)]]j} 

Up 

< Vf >e23 ix,y) = — {[sin[a(x, -y)] - sin[a{-x,y)]]i 

+ [cos[a(a;, -y)] - cos[a{-x,y)]]]} 

< Vf >e24 {x,y) = ^{-[sm[a{x,y)] + sm[a{x, -y)]]i 

+[cos[a{x, y)] - cos[a{x, -y)]]i} {A2) 

with a(x, y) is given by eq. (29). The corresponding relations vahd for other 
regions, can be obteined from (Al) and (A2), using the appropriate rules of 
index and coordinate exchanging (see below). 

Replacing eqs. (Al) and (A2) in (23), we obtain for current density in a 
point of region II , j// 



in{x, y) = [-k(x, y) + \{x, y)]Pi3 + [k(-x, -y) - k(x, y)]Pi2 

-[k(x, y) + l(-x, -y)]P24 + [l(-x, -y) - l(x, y)]P23, (^3) 

where 

k(x, y) = sin a{x, y)i — cos a{x, y)j, 
l(x, y) = sin a{—x, y)\ + cos a{—x, y)j, (^44) 

A cos ^ 

and Pjj = ^^27r ° ^ tanh ( " '^"^ ^ ) . The current density in other regions 
is obtained from j// by applying the following rules of phase and coordinate 
exchanging 

h=hi[{x -y/k,y ^ x/k,k^ 1/k); (i -j,j i); 
iiii=hi[{x y/k,y ^ -x/k,k^ 1/k); (i ^ j,j -i); 

{Pl <P2, P2 <P3, V3 V^4, f4 fl)], (^5) 

jiv =hi[ix -x,y -y,k ^ /c); (i ^ -i, j ^ -j); 

The same relations as (A5) can be used for O^j and < >6»i (the phase 
exchanges have to be replaced by corresponding index exchanges) . 
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Figure Captions 

Figure 1. (a) The superconducting 4-terminal Josephson junction. Four 
coupled superconductind microbridges, going from a week point to the mas- 
sive superconducting banks is the normal resistance of the ith filament 

and ^{T) is the coherence length). 

(b) The mcsoscopic 4-terminal Josephson junction. Four bulk superconduc- 
tors are weakly coupled through a rectangular of two-dimensional electron 
gas (2DEG). 

Figure 2. Dashed line is i — > j trajectory passing through the point p. 
All of this type trajectories are confined in the angle L, W are length 
and width of the rectangular. 

Figure 3. A configuration of the mesoscopic 4-terminal Josephson junc- 
tion. The terminals 1 with 2, and 3 with 4 are short-circuited by the super- 
conducting rings (dashed lines). The phase differences are — (p2 — fi, 'P — 

Figure 4. Absolute value of the induced order parameter | ip{x,y) \ is 
plotted verticaly for the values of phase differences 6' = — V^i = '?r/2, ip = 
cps-cp^^ 37r/2, X = - = 0. The lines of | xpix, y) \= const, are 

shown. 

Figure 5. Vector field plot of the supercurrent density, j(,x,y), inside the 
normal layer. The valuese of phase differences are the same as the Fig. 4. 

Figure 6. The angles Oij for a point in the region //. We have just shown 
^13, ^23, ^24 = ^4 and also the angle a. 
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